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ABSTRACT 

The interaction between Alfven-wave turbulence and the background solar wind affects the cross 
helicity (f cTxv ■ B) in two ways. Non-WKB reflection converts outward-propagating Alfven waves 
into inward-propagating Alfven waves and vice versa, and the turbulence transfers momentum to 
the background flow. When both effects are accounted for, the total cross helicity is conserved. In 
the special case that the background density and flow speed are independent of time, the equations 
of cross-helicity conservation and total-energy conservation can be combined to recover a well-known 
equation derived by Heinemann and Olbert that has been interpreted as a non-WKB generalization of 
wave-action conservation. This latter equation (in contrast to cross-helicity and energy conservation) 
does not hold when the background varies in time. 

Subject headings: solar wind — Sun: corona — turbulence — waves 


1. INTRODUCTION 

Approximately fifty years ago, Parker (1965) and Cole¬ 
man (1968) suggested that waves and turbulence play an 
important role in the heating and acceleration of the so¬ 
lar wind. Since that time, observational, theoretical, and 
numerical studies have produced mounting evidence that 
supports this suggestion. For example, in situ measure¬ 
ments at heliocentric distances exceeding 0.3 AU show 
that turbulent fluctuations pervade the interplanetary 
medium (Goldstein et al. 1995; Bruno & Carbone 2005) 
and that most of this turbulence consists of fluctuations 
that propagate away from the Sun, consistent with a so¬ 
lar origin (Belcher & Davis 1971; Tu & Marsch 1995). 
Remote observations from the Solar Optical Telescope on 
the Hinode satellite reveal the presence of Alfven-wave- 
like motions in the low corona with amplitudes sufficient 
to power the solar wind (De Pontieu et al. 2007). Faraday 
rotation of radio transmissions from the Helios satellite 
are also consistent with theoretical models in which the 
solar wind is powered by an Alfven-wave (AW) energy 
flux (Hollweg et al. 2010). 

In order for AW turbulence to heat the solar wind, AW 
energy that is initially in large-wavelength fluctuations 
must “cascade” to smaller wavelengths, at which the fluc¬ 
tuations can efficiently dissipate. This cascade process 
relies upon the interaction between counter-propagating 
AWs (Iroshnikov 1963; Kraichnan 1965). Because the 
Sun launches only outward-propagating waves, solar- 
wind heating by AW turbulence requires some source of 
inward-propagating AWs. 

One of the most important sources of such inward- 
propagating waves is non-WKB reflection (Heinemann 
& Olbert 1980; Velli 1993; Hollweg & Isenberg 2007). 
Photospheric motions have such long timescales that 
they launch AWs that have radial wavelengths within 
the corona and solar wind that can be comparable to or 
greater than the heliocentric distance. For such waves, 
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the wave phase velocity varies appreciably over one wave 
length, which causes the AWs to undergo partial reflec¬ 
tion as they propagate away from the Sun. 

A number of authors have conducted theoretical and 
numerical investigations of solar-wind turbulence driven 
by non-WKB AW reflection (e.g., Zhou & Matthaeus 
1989; Velli et al. 1989; Matthaeus et al. 1999; Dmitruk 
et al. 2002; Cranmer & van Ballegooijen 2005; Verdini & 
Velli 2007; Chandran & Hollweg 2009; Verdini et al. 2012; 
Perez & Chandran 2013). These authors took the back¬ 
ground solar wind to be steady, and several of them made 
use of a conservation law first obtained by Heinemann & 
Olbert (1980), which Heinemann & Olbert (1980) inter¬ 
preted as a non-WKB generalization of wave-action con¬ 
servation. In this paper, we show that this conservation 
relation does not hold in the case of a time-dependent 
background. We also show that this conservation rela¬ 
tion can be obtained by combining the equations of cross- 
helicity conservation and energy conservation. To the 
best of our knowledge, the equation of cross-helicity con¬ 
servation has not been applied previously to reflection- 
driven AW turbulence in the solar wind. Because the 
equation of “non-WKB wave-action conservation” can 
be obtained from the equation of cross-helicity conserva¬ 
tion, and because cross helicity is conserved regardless of 
whether the background varies in time, the equation of 
cross-helicity conservation in some sense generalizes the 
equation of “non-WKB wave-action conservation” to the 
time-dependent regime. We note that although we allow 
the flow velocity and density to vary in time, our analysis 
is limited to the case in which the background magnetic 
field is fixed. 

In Section 2, we describe and adopt a standard set 
of approximations that has been used in previous treat¬ 
ments of non-WKB reflection of Alfven waves in the 
solar wind. For example, we assume that the back¬ 
ground magnetic held is quasi-radial (Section 2.1) and 
that the fluctuations are transverse and non-compressive 
(Section 2.3). We then summarize the derivation of the 
equations that describe the coupled evolution of the fluc¬ 
tuations and the background flow (Section 2.4). In Sec- 
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tion 3, we describe how the equations in Section 2.4 lead 
to conservation equations for the total energy and total 
cross helicity. We also show how these two conservation 
equations can be combined to recover the aforementioned 
equation of “non-WKB wave-action conservation” when 
the background flow is independent of time. Finally, in 
Section 4, we describe how the equations in our model 
separately conserve the action of inward and outward- 
propagating AWs in the limit of short wavelengths and 
small wave amplitudes. 


2. COUPLED EQUATIONS FOR THE FLUCTUATING 
FIELDS AND BACKGROUND FLOW 

We begin with the equations of ideal magnetohydro¬ 
dynamics (MHD), 
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where p is the mass density, v is the velocity, B is the 
magnetic field, p is the pressure, G is the gravitational 
constant, and Mq is the mass of the Sun. We assume 
that the plasma satisfies the energy equation 
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—pV ■ V — V ■ q, (4) 


where q is the heat flux and 7 is the ratio of specific 
heats. 


2 . 1 . Two Models for the Background Magnetic Field 
We set 

B = Bo + SB, (5) 

where the background magnetic field Bq is a fixed, time- 
independent function of space. We neglect solar rotation 
and consider two models for the background magnetic 
field: 

Model 1: Bq consists of exactly radial magnetic field 
lines filling a region that spans a solid angle of order 
unity as seen from the Sun, with Bq depending only on 
heliocentric distance r. 

Model 2: Bq corresponds to a narrow magnetic flux tube 
in which the magnetic field is nearly but not exactly ra¬ 
dial, which allows for super-radial divergence of the mag¬ 
netic field. In this second case, we take the flux-tube to 
have an approximately square cross section and impose 
periodic boundary conditions on the edges of this cross 
section, as described further in the Appendix. We also 
take the opening angle 9 of the flux tube to be <C 1 . 

In both models, 
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Equation ( 6 ) is exact in model I and correct to leading 
order in 9 in model 2. In both models, we are able to 
discard terms proportional to either V xb or bVb. Such 
terms vanish exactly in the case of model 1 and are much 
smaller than the terms we keep in the case of model 2 . 


2 . 2 . Averages over Surfaces Perpendicular to Bq 

For both background-magnetic-field models discussed 
in Section 2.1, we define the “surface average” of an ar¬ 
bitrary function /, denoted (/), through the equation 

{f) = jJ^dAf, (10) 

where fg dAf denotes an integral of / over a surface S 
that has area A and is everywhere normal to Bq. In 
model 1 , A is the intersection of the modeled region with 
a spherical shell of some radius r. In model 2, the surface 
integral in Equation (10) is confined to the interior of the 
modeled magnetic flux tube. In both models, we are able 
to discard terms of the form (V • aj^) when they arise in 
the derivation of the equations in Section 2.4, where a± 
is everywhere perpendicular to 6 , because such terms are 
much smaller than the terms we keep. We discuss this 
point further in the context of model 2 in the Appendix. 
In model 1, averages of the form (V • aj^) can be dropped 
when a± contains one or more fluctuating quantities in 
part because we assume that 

Lc± < r, ( 11 ) 

where Lc± is the correlation length of the turbulence 
perpendicular to Bq. 


2.3. The Transverse, Non-Compressive Approximation 

As in Equation (5), we set each of p, v, and p equal 
to the sum of a background value (denoted with a ’ 0 ’ 
subscript) and a fluctuating part {5p, 5v, and 5p). We 
define 

V = ti||5 + ^v_L, (12) 


where -6 = 0. Eor the scalar quantities U||, p, and 
p, we define the background quantities to be surface av¬ 
erages: 
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and 
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We refer to Equations (14) through (17) as the trans¬ 
verse, non-compressive approximation. Observations 
provide some support for this approximation. For ex¬ 
ample, in situ measurements show that turbulent fluc¬ 
tuations in the solar wind are weakly compressive and 
preferentially transverse (see, e.g., Tu & Marsch 1995; 
Goldstein et al. 1995; Bruno & Carbone 2005). Radio¬ 
scintillation observations further indicate that (5/9 <C po at 
heliocentric distances as small as a few solar radii (Coles 
& Harmon 1989; Markovskii & Hollweg 2002; Chan- 
dran et al. 2009). On the other hand, transverse, non- 
compressive fluctuations nonlinearly generate compres¬ 
sive fluctuations and longitudinal fluctuations at some 
level, an effect that we neglect. Thus, while the trans¬ 
verse, non-compressive approximation may apply to the 
bulk of the fluctuation energy, the equations we derive 
in Section 2.4 do not account for all of the physical pro¬ 
cesses occurring in solar-wind turbulence. 


Averaging Equation (1) as in Section 2.2, we obtain 
^ = -V • (poC/S). (22) 


Subtracting Equation (22) from Equation (1), we find 
that 
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Because of the transverse, non-compressive approxima¬ 
tion (Equations (14) through (17)), the density fluctua¬ 
tions have no effect on the flow to leading order. Equa¬ 
tion (23) thus describes the evolution of passive-scalar 
density fluctuations in the expanding solar wind. 

We define a normalized magnetic fluctuation 
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2.4. Coupled Equations for the Fluctuating and 
Background Quantities 

In this section, we present the equations that describe 
the fluctuations and background flow to leading order 
in 0, Lcj_lr, Spfpo, and 6p/po. To obtain an equation 
describing the average parallel velocity, we take the dot 
product of Equation (2) with b and then average the 
resulting equation over a surface perpendicular to b as 
described in Section 2.2. This yields 


and the Elsasser variables 

T (25) 

Given our sign convention in Equation (25), z'^ iz~ ) rep¬ 
resents non-compressive, Alfven-wave-like fluctuations 
that propagate in the direction of Bq {—Bq). By com¬ 
bining Equations (19) and (21), we find that 
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The term proportional to a on the right-hand side of 
Equation (18) is the radial component of the averaged 
MHD Reynolds stress {SB ■ VSB/iir — poSvj_ ■ VSvj_), 
and was obtained previously by Usmanov et al. (2011, 
2014). 

Upon multiplying Equation (18) by b and subtracting 
the resulting equation from Equation (2), we obtain the 
perpendicular momentum equation, 

Po§-^Sv± + poU ■ VSvA + ^^SvA = -Va {SB) 


where 


Bo 


(27) 


is the Alfven velocity. Equation (26) was previously used 
by Chandran Sz Hollweg (2009) and is a specialized form 
of the more general Elasser-variable equation obtained 
by a number of authors (e.g., Zhou & Matthaeus 1990; 
Velli 1993; Verdini & Velli 2007; Zank et al. 2012), in 
which we have used Equation (6) to replace the quantity 
z^ ■ V{—U ± va) appearing in those studies with the 
quantity {af2)z^{—U ± va)- We rewrite Equation (26) 
in terms of the Elsasser stream functions and Elsasser 
vorticities in the Appendix. 
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3. CONSERVATION LAWS 

To obtain an equation expressing conservation of total 
energy, we first take the dot product of Equation (26) 
with 2z^ and average over a surface perpendicular to b 
to find 


(5n is the fluctuating part of p -I- B'^/Stt, the quantity 
(a)j_ is defined via Equation (7) for arbitrary a, and 
Vj_./ = (V/)i for an arbitrary scalar function /. The 
projection of Equation (3) onto a plane perpendicular 
to b yields 


|((0±)2) + (J7 ± ua)-V((z±)2) = -aiUT va) {z+-z-) 
±{{z^)^-z+-z-)(v-VATl^-Uy (28) 


= Bo ■ VSv± - U ■ VSB + ^{USB - BoSv±) 
-SBV ■ U + {SB ■ VSv± - • VSB)^ . 


We then take the sum of the following equations: Equa¬ 
tion (18) multiplied by U; Equation (22) multiplied 
by C/^/2; the “plus version” of Equation (28) multiplied 
by po/4; the “minus version” of Equation (28) multiplied 
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by /9 o/ 4; and the average of Equation (4) over a surface 
perpendicular to b. This yields 
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is the energy density of the turbulent fluctuations, 
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is the gravitational potential, and 
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is the surface-averaged total-energy flux. 

The surface-averaged cross-helicity density is 


(33) 
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We obtain the equation expressing total-cross-helicity 
conservation by adding the following equations: the “mi¬ 
nus version” of Equation (28) multiplied by y^7f/io/2; the 
“plus version” of Equation (28) multiplied by —y/'!Tpof2; 
and Equation (18) multiplied by Bq/po- This yields 
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is the surface-averaged cross-helicity flux, and 

h{r)=j'dr—^. (37) 

J Po dr 

The fact that cross helicity is conserved in the pres¬ 
ence of non-WKB wave reflection is perhaps surprising. 
Most studies of incompressible MHD turbulence focus on 
the case of a stationary background, in which the cross 
helicity arises entirely from the turbulent fluctuations. 
In that case, converting z^ fluctuations into z^ fluctu¬ 
ations would violate cross-helicity conservation. In con¬ 
trast, the interaction between fluctuations and a moving, 
inhomogeneous, and time-dependent solar wind changes 
the cross helicity via two mechanisms: non-WKB wave 
reflection, which changes the cross helicity in the fluctu¬ 
ations, and the transfer of momentum from the fluctua¬ 
tions to the background plasma, which alters the cross- 
helicity content of the background flow. The combined 
effect of these two mechanisms conserves the total cross 
helicity in the sense of Equation (35), which, when in¬ 
tegrated over some arbitrary volume, implies that the 


change in the total cross helicity within that volume 
equals the amount of cross helicity that flows into that 
volume through its boundaries. 

We note that total cross helicity is also conserved in 
weak, homogeneous, compressible MHD turbulence, de¬ 
spite the fact that interactions between Alfven waves 
and magnetosonic waves convert z^ energy into z^ en¬ 
ergy (Chandran 2008). In that problem, there is no flow 
of cross helicity through the boundaries, and the change 
in the cross helicity of the fluctuations is exactly offset by 
the change in the cross helicity of the background. The 
cross helicity of the background changes because the reso¬ 
nant three-wave interactions that convert z^ fluctuations 
into z^ fluctuations simultaneously generate a small, av¬ 
erage, background flow parallel or anti-parallel to Bq. 

We can combine the equations of cross-helicity con¬ 
servation and energy conservation by first multiplying 
Equation (35) by pqU/Bq and then using Equation (29) 
to rewrite the term {poU/B)V■ (Boh) = U-Vpo in terms 
of variables other than pq. Recalling that Bq and U are 
parallel, which implies that Bq ■ V{poU/Bo) = —dpo/dt, 
we find after some algebra that 
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When the background plasma is steady, the right-hand 
side of Equation (38) vanishes, and Equation (38) is 
equivalent to Equation (26) of Heinemann & Olbert 
(1980), which those authors interpreted as the gener¬ 
alization of AW action conservation to the non-WKB 
regime. Although Heinemann & Olbert (1980) derived 
their Equation (26) for linear waves, their Equation (26) 
is also valid in the nonlinear regime, provided dpo/dt = 
dU/dt = 0, as can be seen from Equation (38) above. 
On the other hand, when the background plasma varies 
in time, the right-hand side of Equation (38) is in gen¬ 
eral nonzero. Thus, Heinemann & Olbert’s (1980) Equa¬ 
tion (26) does not extend to the case of a time-dependent 
background. 


4. WAVE ACTION 

Bretherton & Garrett (1968) considered the propaga¬ 
tion of linear waves in slowly varying, inhomogeneous, 
moving media in the WKB limit of short wavelengths 
and short wave periods. They took the waves to satisfy 
a dispersion relation of the form 

w = D(fe, A(r,t)), (39) 

where A(r,t) is some slowly varying function of position 
and time. The group velocity of the waves is then 

c = VfeD, (40) 

where Vfc denotes the gradient operator in wavenumber 
space, and the frequency varies along a ray path accord¬ 
ing to the equation 

duj dn dX 
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Bretherton & Garrett (1968) showed that for a wide 
class of conservative systems, including Alfven waves in 
a time-dependent, inhomogeneous, moving medium, 
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+ (V • c) 
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where c is the group velocity, 

oj' = oj-k-U (44) 


is the wave frequency measured in the local rest frame of 
the medium, U is the velocity of the medium, and fw is 
the energy density of the waves. 

For the case of WKB Alfven waves propagating in a 
radial magnetic field, £w = po(( 2=‘=)^)/4, n = kr{U ± 
va), a = t/ ± VA, and uj' = ±krVA- Upon multiplying 
Equation (43) by uj and making use of Equation (41), we 
obtain 


Jt v~ ) 
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To see how Equation (45) is recovered as a limit¬ 
ing case in our analysis, we multiply Equation (28) by 
Po{U ±va)/(4ua) and simplify the resulting expression 
using Equation (22) and the identities Bo-V{poU/Bo) = 
-dpofdt and V • {pqVa) = -poV • va = (1/2)ua ■ Vpo- 
After some algebra, we obtain 
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In the limit of short wavelengths and small wave am¬ 
plitudes, ( 2 :+ • 2 ;“) —>• 0. In this limit, Equation (46) 
reduces to Equation (45). When dpo/dt, dU/dt, and 
dvA/dt vanish, subtracting the “minus version” of Equa¬ 
tion (46) from the “plus version” of Equation (46) repro¬ 
duces Equation (38) with the right-hand side of Equa¬ 
tion (38) replaced by zero. 


5. CONCLUSION 

Conservation laws play a fundamental role in the study 
of turbulence, because they are among the few analytic 
results that can be used to gain insight into the physics 
of turbulent systems. Eor example, energy conservation 
underpins the concept of an energy cascade, in which 


nonlinear interactions among fluctuations transfer fluc¬ 
tuation energy in a loss-free manner from large scales to 
small scales. This idea ultimately explains why the tur¬ 
bulent heating rate can be determined solely from the 
properties of the turbulence at large scales (i.e., the in¬ 
ertial range or the outer scale), regardless of the mecha¬ 
nisms that dissipate the energy at small scales. Conserva¬ 
tion of magnetic helicity in MHD turbulence leads to the 
concept of an inverse cascade of magnetic helicity, which 
plays an important role in turbulent dynamos (Erisch 
et al. 1975; Pouquet et al. 1976). In this paper, we have 
shown that a third conservation law, that of cross helic¬ 
ity, applies to non-WKB AWs and reflection-driven AW 
turbulence in the solar wind. This result is in some ways 
surprising, because non-WKB reflection converts 2 ^ fluc¬ 
tuations into fluctuations, thereby altering the cross- 
helicity content of the fluctuations. The total cross he¬ 
licity is nevertheless conserved because the fluctuations 
exert a force on the background solar wind, which alters 
the cross-helicity content of the background flow. 

Our finding that cross helicity is conserved by non- 
WKB AWs and reflection-driven AW turbulence is im¬ 
portant for a few reasons. Eirst, it implies that cross 
helicity can be exchanged between the fluctuations and 
the background flow without loss. Second, in contrast 
to the equation of “non-WKB wave-action conservation” 
derived by Heinemann & Olbert (1980), cross-helicity 
conservation holds even when the background flow varies 
in time. (Also, as discussed in Section 3, the equations of 
cross-helicity and energy conservation can be combined 
to recover Heinemann & Olbert’s (1980) conservation law 
when the background solar wind is time-independent.) 
Third, the coupled equations for the fluctuations and 
background flow in Section 2 can be solved numerically 
to provide new insights into the heating and accelera¬ 
tion of the solar wind by reflection-driven AW turbu¬ 
lence, and cross-helicity conservation provides a valu¬ 
able benchmarking tool for such simulations. There is 
a growing interest in numerical simulations of the so¬ 
lar wind that incorporate AW turbulence (see, e.g., Us¬ 
manov et al. 2011, 2014; van der Holst et al. 2014), in 
part because of the upcoming launch of Solar Probe Plus. 
This pioneering mission will shed new light on the mech¬ 
anisms that heat and accelerate the solar wind by pro¬ 
viding the first-ever in-situ measurements of the solar- 
wind acceleration region. By providing new insights into 
one such mechanism (reflection-driven AW turbulence) 
as well as a valuable benchmarking tool for certain types 
of numerical codes, our results may ultimately contribute 
to a deeper understanding of the solar wind’s origin. 
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APPENDIX 

THE NARROW-MAGNETIC-FLUX-TUBE APPROXIMATION 

As discussed in Section 2.1, our results apply under either of two different assumptions about the geometry of the 
background magnetic Held. The background field can be either exactly radial throughout a region spanning a solid 
angle of order unity as seen from the Sun, or the background field can be approximately radial within a narrow magnetic 
flux tube centered on a radial magnetic field line. In this appendix, we consider the case of a narrow magnetic flux 
tube. We work in spherical polar coordinates and take 0 = 0 to correspond to the exactly radial background magnetic 



6 


field line that coincides with the axis of the magnetic flux tube. We restrict our analysis to a region within which 

0<1. (Al) 

We assume that 

(A2) 


dB^r dBoe „ 
= ... ■ = ... ■ = U 


d(j) d(j) 


and define 


H{r) = \Bo{r,e = 0)\. (A3) 

The condition that V ■ B = 0 implies that dBg/dO = 0 at 0 = 0. We require that V^Bq be finite, which implies that 
dB[)rld9 vanishes at 0 = 0, so that 

B^r(.r,e) = H{r)[l + 0{e^)\. (A4) 

The condition V • Bq = 0 then implies that 


^00 


We assume that 


Equations (A2), (A4) and (A5) imply that 
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2r dr 
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r dH 
H dr 
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r - eJ^^ir^H) I [1 + 0(02)] , 


2rH dr 

where b = Bq/Bq. It follows from Equation (A7) that 

V • 6 = cr + O 
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where 

It also follows that for any vector a, 

a ■ Vb 


1 dH 
H dr ’ 


|a, + 0 
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(A7) 

(AS) 
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(AlO) 


where a|| = b ■ a and aj_ = a — a||b. 

The quantities V x 6 and b ■ Vb relate to twist and curvature of magnetic field lines and satisfy the relations 


and 


where 


b-Vb^V^-tp + O 



Vx S = bx (S-VS), 


V' 


02 r 


2 —(or) + (or 
dr 



(All) 

(A12) 

(A13) 


Equation (A12) is not a vector identity, but is exact because of Equation (A2). It follows from Equations (A6), (A8), 
(All) and (A12) that 

|V X S| ~ |S-VS| ^ O < 

At several points in Sections 2 and 3, we average equations over surfaces that are everywhere perpendicular to JBq, 
as described in Section 2.2. To specify these surfaces mathematically in the case of a narrow magnetic flux tube in 
which Bq is not exactly radial, we introduce the vector potential Aq associated with the background magnetic field 
and define Clebsch coordinates (Euler potentials) a and /3 that are related to Aq through the equation Aq = aV/3, 
which yields 

Bq = Vax V/3. (A15) 


|V- b| - O 


(AM) 



7 



Fig. 1.— Clebsch coordinates (o:,/5,s). Surfaces of constant s are perpendicular to Bq. The values of o: and (3 are constant along field 
lines of Bq. 

Since Bq ■ Va = Bq ■ V/3 = 0, a and f3 are constant along the magnetic field lines of Bq. The particular Clebsch 
coordinates that we use are 

a = [H{r)Y/^x (A16) 


l3 = (A17) 

where {x,y,z) are Cartesian coordinates, and the positive 2 axis coincides with 9 = 0. When Equations (A16) and 
(A17) are substituted into Equation (A15), the resulting value of Bq satisfies Equations (A4) and (A5), as required. 
We introduce a third coordinate s such that surfaces of constant s are perpendicular to Bq, with s = r at 0 = 0: 




The {a,(3,s) coordinate system is illustrated in Figure 1. 

For any function f(a,j3,s), we define (/) to be the average of / over a surface of constant s with —5 < a < 6 
and —S < (5 < 5, where d is a constant that fixes the width of the flux tube. Two such surfaces are sketched with 
cross-hatched lines in Figure 1. Thus, 


where 


A = da dp — 

J-S J-S Bq 


is the area of the averaging surface. 

Equations (A4), (A5), and (A18) imply that |i3o| = H{s)[l -I- 0{9‘^)]. We can thus rewrite Equation (A19) as 


^ d/3/[I+ 0(0^)] 


where 


$ = 4d^ 


is the magnetic flux through the averaging surface, which is independent of s. Equation (A18) implies that |Vs| = 
1 -I- 0{9‘^), so that 


It follows from Equation (A21) that 


S-v/ = ^[1 + 0(02)]. 





where 9o is the value of 9 at the middle of one of the edges of the averaging surface - e.g., at coordinates (a,/3, s) = 
(<5,0,s). 

We say that a scalar function g satisfies periodic boundary conditions in the plane perpendicular to Bq if it obeys 
the relations 


g{6,l3,s) = g{-6,l3,s) g{a,6,s) = g{a,-S,s) 


(A25) 


for all a £ (—S,S), all /3 G (—d, d), and all s. We say that a vector a satisfies periodic boundary conditions in the 
plane perpendicular to Bq if the vector’s components a ■ Va, a ■ V/3, and a ■ Vs satisfy the periodicity relations in 
Equations (A25). If a vector a± is periodic in the plane perpendicular to Bg and satisfies a± ■ Bg = 0 everywhere, 
then Stokes’ theorem can be used to show that 


(V • a^} 


r 


(A26) 


where we have assumed that the characteristic length scale of aj_ perpendicular to Bg is ~ 9gr. The rms value of 
V-a± on the averaging surface is ~ (|aj_p)^/^/(0o?')- Equation (A26) thus implies that the average of V-aj_ is reduced 
relative to its rms value by a factor of ~ 9q. (This reduction factor would be even smaller if the characteristic length 
scale of a_L perpendicular to Bg were much smaller than 9gr.) This reduction enables us to drop averaged quantities 
of the form (V • a±) in Section 2, because they contribute only higher-order corrections to the equations presented. 


ELSASSER STREAM FUNCTIONS AND VORTICITIES 
We define the Elsasser stream functions and the Elsasser vorticities through the equations 


= bx VC 


= vic^, 


where 


Vi/ = V • (Vi/) = V • V/ - S(S • V/) 


(Bl) 

(B2) 


for any function /. By taking the cross product of Equation (26) with b and then taking the divergence of the resulting 
equation, we can rewrite Equation (26) in the form 

^r!± + ([/±UA)(h-VfI± + |fI±) = |(-C/±z;A)fI^ + i (^V-iiaT^V-C/^ (fl+-fl-) 

({r, vic+} + {c+, viri ± vi{c-,c+}), (bs) 

where 

{/, 5 } = 6.(Vi/x Viff). (B4) 

Here, we have assumed that either 0 <C 1 for the case in which the background magnetic field corresponds to a narrow 
magnetic flux tube or Lc± <C r for the case in which Bg is radial throughout a region spanning a solid angle of 
order unity, where Lc_l is the correlation length of the fluctuations perpendicular to Bg. Equation (B3) generalizes 
Equation (A4) of van Ballegooijen et al. (2011) to account for the background flow U. Equation (B3) is the same 
equation that was solved numerically by Perez & Chandran (2013). (Note that the minus sign on the right-hand side 
of their Equation (10), which was erroneous, was a typo in their paper, not an error in their code.) The form of the 
nonlinear term on the last line of Equation (B3) is the same as in the RMHD equations derived by Schekochihin et al. 
(2009), except that our corresponds to their In the homogeneous-background limit (in which a, V ■ va, and 
V • U vanish). Equation (B3) reduces to Equation (21) of Schekochihin et al. (2009). 
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